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We consider an integrability test for ultradiscrete equations based on the singularity confinement 
analysis for discrete equations. We show how singularity pattern of the test is transformed into that 
of ultradiscrete equation. The ultradiscrete solution pattern can be interpreted as a perturbed solu- 
tion. We can also check an integrability of a given equation by a perturbation growth of a solution, 
namely Lyapunov exponent. Therefore, singularity confinement test and Lyapunov exponent are 
related each other in ultradiscrete equations and we propose an integrability test from this point of 
view. 



I. INTRODUCTION 

Integrability is an important concept in nonlinear equations. If a given equation turns out to be integrable, we can 
get many exact structures from the system, for example, conserved quantities, symmetries, exact solutions, and so on. 
Therefore, it has been an important problem to test integrability of equations. For differential equations, it is well 
known that the Painleve test is powerful to detect integrability. [I]] However, for difference equations, the Painleve test 
cannot be applied due to discreteness of independent variables, and the singularity confinement (SC) test is proposed 
instead. Q 

Let us consider the following multiplicative type of difference equation 

Xn+xX^Xn-l = a\ n X n + 1. (1) 

This equation for a = 0, 1, 2 with A = 1 belongs to so-called Quispel- Roberts-Thompson (QRT) system ||Q which is 
a large family of integrable second order ordinary difference equations. The term "integrability" is somewhat more 
delicate in discrete system than continuous one. For QRT system, it is integrable in a sense that it has a quartic 
conserved quantity and thus the general solution is expressed by elliptic functions. For generic A, equation (Q) with 
a = 0, 1,2 are known as discrete analogues to Painleve I equation, which are considered to be integrable in a sense 
that it passes the SC test. || It is considered to be non-integrable for other a. 

The SC test is applied to equation (P as follows. Let us assume a = 2 and A = 1 to make discussions easier. If 
initial x's(assume xo and x\) are xq = / (non-zero finite) and x\ — — — respectively, then we have X2 — 0, and X3 
becomes singular. In order to see the behaviour of this singularity, we introduce a small parameter 8 (~ 0) and put 
xq = f and x\ = — — + 5. Then, successive iteration and asymptotic evaluation in S gives a singularity pattern as 
follows: 
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The above pattern shows; (a) A singularity due to X\ occurs at n = 3, (b) The singularity is confined, that is, does 
not spread on the whole lattice, (c) Information on xq pass through the singularity to x n at n > 4. Following to SC 
test, both locality of singularity and preservation of information on initial data strongly indicate an integrability of 
the equation. 

There is another singularity pattern using initial data xq = f and x\ =6: 
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The above SC pattern shows; (a) xg n = finite and x^n+i — 0(8) for any n, (b) Information of £0 pass through 
periodic singularities to x n at large n. In this pattern, singularity spreads on the whole lattice but singular points 
are confined by finite values which have information on initial Xq . Although such periodic singularity pattern has not 
been discussed anywhere to the authors' knowledge, we may consider that this pattern also suggests integrability. 
Next let us consider an ultradiscrete analogue to equation (|l]). P pp|Jio|-[T4| Applying the following transformation, 
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(4) 



we obtain 

X n+1 +aX n + X n - 1 = £log(l + e ^" + " i+A)/e ). (5) 
If we take a limit e — > +0, we get 

X n+1 +aX n + X n -i = max(0, X n + n L + A), (6) 

using a formula 

lim e log{e a/e + e b/e + ••■)= max(a, b, ■ ■ ■). (7) 

If initial data Xq, X\ and parameters L and A of (^) are all integer, Jf„ for any n is always integer. In this sense, the 
dependent variable of equation (||) is discretized through a limit (0). The above discretization process on dependent 
variable is called 'ultradiscretization' and equation (|^) is an ultradiscrete analogue to equation. (|l|) M We obtain a 
cellular automaton from an ultradiscrete equation, if it is possible to restrict values of dependent variable in a range 
of integers. Thus, it is expected that ultradiscrete equations connect the continuous and the discrete worlds, that is, 
differential or difference equations and cellular automata. 

A large difference between ultradiscrete equations and discrete ones is an existence of singularities. Since ultradis- 
crete equations are piecewise linear, singularities in normal sense do not exist. Therefore, we cannot use SC test as it 
is. 

In this letter, we demonstrate our trial to detect integrability of ultradiscrete equations based on the SC test for 
discrete equations. In the next section, we show how SC test is transformed in equation (^). Then, we discuss 
differences between ultradiscrete analogous to integrable difference equations and those to nonintegrable ones. In 
section 3, we show that SC test is related to Lyapunov exponent in ultradiscrete equations. Lyapunov exponent is a 
growth rate of perturbation of solution, by which we can test integrability of a given equation. Jig ] In the last section, 
we give concluding remarks. 



II. ULTRADISCRETE ANALOGUE TO SINGULARITY PATTERN 

In this section, we discuss how the singularity pattern of the difference equation (Q) is transformed to that of 
ultradiscrete analogue (|J). First we note that the case of a = 0,1,2 in autonomous case (L = 0) of equation (||) 
should be considered to be integrable, since; (a) Each of them has a conserved quantity (b) Each of them admits 
a general solution which is expressed by ultradiscrete analogue of elliptic functions. JTo[ | Therefore, we can examine 
whether a transformed singularity pattern is useful to test an integrability of ultradiscrete equation or not by comparing 
the pattern for a = 0, 1, 2 with that for other cases. 

When we ultradiscretize equation ([!]) to obtain equation (^J), we use transformation of variables (Q). In this 
transformation, if L and A are not 0, this means that we introduce singularities into parameters a and A as well as x 
in equation ([[]). In order to discuss effects of singular parameters separately, first we consider a case of finite a and 
A, that is, A = L = 0. Moreover, let us assume a = 2 case, 

X n+ i + 2X n + X n -i = max(0, X n ), (8) 

as a concrete example. In the transformation (Q), x (also a and A) must be always positive. Therefore, we cannot 
obtain an ultradiscrete pattern corresponding to pattern (||) because x n necessarily becomes negative at certain n's. 

On the other hand, if / and 5 are positive, x n in the singularity pattern (|3|) is always positive and can be ultradis- 
cretized. In the limiting process, relation between those values and ultradiscretization parameter e is important. Let 
us assume that / is finite independent of e and <5 = e~ K l e where K is positive and finite. Then, we get the following 
pattern on X n from pattern (]|) through ([|). 
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This pattern is strictly periodic with period 8. Note that this pattern can also be obtained from equation (|8|) with 
initial data Xq — and X± = —K directly. We can easily see that positive, zero, negative X n correspond to infinite, 
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finite, infinitesimal x n respectively. Periodic confinement of singularities of x n is transformed into separation of 
positive values by O's in X n . However, only K corresponding to 8 survives and information on / is lost in the pattern. 

It is important in the SC test whether the information of initial data survives or not. Therefore, next we take 
/ = e p / £ and 8 = e~ K ^ £ where p is finite and \p\ <C K. The initial data become Xq — p and X\ = —K, and the 
following pattern is obtained from pattern (^), 
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This pattern is also periodic with period 8. The information of initial data, namely p, clearly survives for large n. 
It is easy to see that we can obtain this pattern directly by successive iteration of equation (g). Therefore, we can 
conclude that the equation and its singular solution pattern are consistently transformed from difference equation to 
ultradiscrete one at least in this example. Similar results are obtained for the other integrable cases (a = and 1). 

Next, we test non-integrable case of equation (Q), a = 3. Let us assume A = 1 for simplicity, then equation (|l|) 
becomes 



Xn+\X n X n —\ 
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and we get a singularity pattern 
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with xo = f and x\ = S. We can easily prove that x n = 0(8 Pn ) and p 2n ~ —^((2 + \/3)" — (2 — V3) n ) and 

P2n+i ~ (2 + \/3) n — ^2~ X (2 — \fi) n . Singularities are not confined and information on initial data is lost for 
large n when 5^0. Following to SC test, this strongly implies that equation (jll]) is not integrable. 
The corresponding ultradiscrete equation is 



X n+1 + 3X n + X n -i = max(0, X n ). 



(13) 



and a pattern becomes 
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In this pattern, the first term including K in every X n (n > 1) becomes dominant and information on p becomes 
negligible. 

In the above examples, we have shown singularity pattern of difference equation (Q) for finite a with A = 1 
and corresponding pattern of ultradiscrete equation (H) simultaneously. The values / (finite) and 1/8 (singular) 
in equation (|l|) corresponds to p and K (|p| <C K) in equation (||), respectively. When equation (0) is integrable 
(er = 0, 1,2), information of p survives for large n in corresponding ultradiscrete equations. Confined singularity in 
the SC test on difference equation ([!]) corresponds to separation of positive values including K by small p. 

On the other hand, when equation (0) is non-integrable (er = 3), information on p is negligible and K becomes 
dominant for large n in corresponding ultradiscrete equation, and separation of positive values by p does not occur. 

We note that the above singular solution patterns for ultradiscrete equation are obtained without inconsistency. 
Namely, the pattern obtained by taking ultradiscrete limit on the singularity pattern of difference equation always 
coincides with that obtained by successive iteration of ultradiscrete equation directly. 

Therefore, we could obtain an interpretation of the SC test on the ultradisctete equation (^) with L = A = 0, and 
observe a critical difference between the case of a — 0, 1, 2 and other cases. 
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III. LYAPUNOV EXPONENT OF PERTURBED ULTRADISCRETE SOLUTION 



So far, we have considered the ultradiscrete equation (|^) with L = A = 0. However, the parameters L and A can 
effect a solution drastically. For example, let us take a — 2 case in equation (0) again, 

X n+1 +2X n + X„_i = max(0, X n + n L + A), (f 5) 

and put L = and A = 1. Let us start with the initial data Xq — and X% = — 1 . Then we see that this solution 
is periodic with period 3. However, starting with Xq = and X\ = —2, we see that the period is now 20. Compare 
this phenomenon with a fact that the period for L = A = was 8 for any Xq = p and X\ = —K. Therefore, the 
period of solution can change by the combination of parameter A and initial data Xq, X\. Moreover, in a case of 
L ^ 0, equation ( |l5|) becomes non-autonomous equation. Therefore, since the values of X n 's grow due to the term 
nL, both criterions obtained in the previous section do not work as they are. That is, information on initial value 
becomes negligible and separation of positive values by p does not occur. A solution to equation (|l5| ) in the case of 
A = 3, L = 2, X = 0, Xi = 1 is shown in Fig. 0. 

These observations imply that we cannot simply transform a singularity pattern of difference equation to ultra- 
discrete one for L ^ and 4^0. This is because taking the parameters L and A to be non-zero means that we 
have introduced singularities on parameter A and a in the difference equation ([l]), respectively. Therefore, we must 
improve an integrability test for ultradiscrete equation. 

Let us look back the singular solution patterns ( ^p| ) and ( |l4| ) of the ultradiscrete equations (0) and ( |l3| ) , respectively. 
We can give another interpretation of p and K . if is a parameter giving a solution orbit with initial data Xq = 0, 
X\ = —K and p is a perturbation to the orbit. From this viewpoint, growth rates of perturbation in the patterns 
( |To| ) and (|l4| ) are clearly different. Therefore, it might be possible to regard SC analysis on ultradiscrete equations is 
nothing but analysis on a growth rate of perturbation. 

Let us assume that X n is a solution to equation ( [l5| ) and X' n is one perturbed by p. Then, amplitude of perturbation 
is defined by a n where a n = \ "~ " |. Figure |2| shows a n for A — 3, L = 2, Xq = 0, X = p and X\ — X[ — 1. We 
can see that a n grows linearly for n < and does not grow for n > 0. Therefore, the solution is not chaotic and we 
can easily estimate a global behavior of the solution. 

A quantitative index to check integrability of a system is Lyapunov exponent which is a mean growth rate of 
perturbation. [jl5| If we define A„ by ^loga n , Xoo (A_oo) is the exponent. When Aoo > 0, the perturbation grows 
exponentially and the system becomes chaotic. As for equation ([13|), Aio n for Xq = 0, X' = p and X\ — X\ = 1 
is plotted in Fig. |[ A„ is asymptotically about 0.66 for large n. Therefore, we can conclude that equation (|l3|) is 
chaotic and this result is compatible with that of the previous section. 

Next we try another example, the cutoff Toda equation [[l6| 

x n+l = x n/y n 

y n +i =y n (h + x n+1 )/(h + l/x n+1 ) ' v ; 

where a and h are constant parameters. Eliminating y, this equation reduces to the following ordinary difference 
equation on x n , 

x n+ ix n -i =x„~ <T (l + hx n y /(h + x n y . (17) 



Equation (17) with a — 1 and that with a = 2 are nothing but the autonomous discrete Painleve H and HI equations, 
respectively. [p| According to SC test, only the cases of a = 1 and 2 are integrable and not otherwise. Here we discuss 
details of perturbed solutions of ultradiscrete analogue to ( |l6| ) in the cases of a = 2 (integrable) and 3 (non- integrable) . 
If we introduce x n = exp(A„/e), h — exp(H/e) and take the limit e +0, we obtain 

X n+1 + X„_i = (2 - a) X n + g (max(0, X n + H) - max(F, X n )), (18) 

from equation (jlj]). If H ^ 0, this means that equation ( |Tt| ) includes singular parameter h. In the case of a — 2, we 
see that the difference equation (|l7|) has the following conserved quantity, 



x n+1 x n + 2h(x n+1 + x n ) + /i 2 (^±i + -^-) + 2h{-^— + —) + . (19) 

Xn Xfi-\-\ X n X n -\-\X n 

Then, we obtain a conserved quantity of ultradiscrete equation ( p"8| ) with a = 2, 

max(A„ + X n+ i, —X n — X n+ i, 2H — X n + X n+ i, 2H + X n — X n+ i), (20) 
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through the above ultradiscrete limit. Define K by max(Xo + Xi, —Xq—Xi, 2H — Xq + X%, 2H + Xq — X\), then any 
point (X n , X n+ i) (n > 0) in phase space always exists on sides of a rectangle with vertices (H — K, —H), (H, K — H), 
(-H, H — K), (K — H, H). If we use perturbed initial data {Xq + pi, X\ + P2) where p\ ~ and pi ~ 0, then an 
orbit of (X n , X n+ i) shifts infinitesimally but it is still stable. 

To estimate a growth of perturbation, let us calculate amplitude of perturbation a n by solutions X n and X' n which 
are the solutions from given initial data (Xq, Xi) and (Xq + p, Xi), respectively. If we choose H — 5, Xo — and 
X\ = 1, then X n becomes periodic with period 24 and <Z24 n = 20n — 1. If we take H — 10, Xq = 11 and X\ = 17, 
period is 9 and ag n = 5n + 1. In both cases, a n grows linearly on n. These observations suggest that equation Jl8| ) 
with a = 2 is integrable. 

On the other hand, in the case of a = 3, the growth of perturbation changes drastically. Choosing H = 5, Xo = 
and Xi = 1, X n becomes periodic with period 18 and a Wn = 1, 901, 326101, 118047661, 42732927181, • • •. Clearly a n 
shows an exponential growth and does not grow linearly on n. If we take H — 10, Xq — 11 and X\ = 17, period is now 
132 and ai 32n = 1, 84020277977, 6297688959151021350515, 472039455011914887138396062962253, • • •. Fig. | shows 
a figure of the solution orbit and Fig. |5| shows a relation between n and log 0132,1, both in the latter case. From Fig. ||, 
we can see that the orbit is more complicated than that in a = 2 case. From Fig. ||, ai32n ~ e 25n and A ~ 25/132. 
Therefore, we can consider that solutions to equation ( |l8| ) with a = 3 become chaotic and thus equation (|l^) is not 
integrable. 

IV. CONCLUDING REMARKS 

In this letter, we obtained the following results. 

(a) If parameters of ultradiscrete equation are all 0, that is, if corresponding parameters of difference equation are 
not singular, we can consistently transform singularity pattern of difference equation into ultradiscrete one. 
Finite value and singularity of solution to difference equation correspond to perturbation and finite value of 
ultradiscrete solution respectively. From other point of view, the SC analysis on difference equation can be 
regarded as perturbation analysis of solution on ultradiscrete equation. 

(b) If parameters of ultradiscrete equation are not 0, that is, if corresponding parameters of difference equation 
are singular, we cannot consistently transform singularity pattern of difference equation into ultradiscrete one. 
However, we can observe a growth rate of perturbation, that is, Lyapunov exponent and check integrability by 
the exponent. 

Moreover, we point out an advantage of the above test. In SC test on difference equation, we need asymptotic 
evaluation of a solution including small parameters. It often needs a large amount of symbolic manipulation and 
becomes hard even with help of a computer. If we can transform such a difference equation into ultradiscrete one, we 
easily calculate a solution pattern. Therefore, when we test integrability of difference equation, it is much easier to 
test a corresponding ultradiscrete equation instead. 

Finally, we give future problems. 

(i) Since integrability test described in (b) is generic, we can use it to any ultradiscrete equation. It is easy to 
conclude a given equation is chaotic if Lyapunov exponent of a solution from a particular initial data is positive. 
However, it is difficult to judge a given equation is integrable. Because we must check Lyapunov exponent is 
not positive for ANY initial data. To avoid this, we need a rule to select initial data to make the test finite. 

(ii) As mentioned in the introduction, we get some cellular automata using ultradiscretization when we can restrict 
values of dependent variables in a range of integers. []l3| We can test integrability of the cellular automata 
using perturbation. Therefore, intermediate values as well as discrete values are significant. However, cellular 
automata are normally defined by binary operations irrelevant to difference equations and we cannot obtain 
information on intermediate values. Thus, further considerations are required to study cellular automata using 
our approach. 

(iii) It is remarkable that Hietarinta and Viallet reported that there are some difference equations that passes the 
SC test but not integrable. E^] Therefore, it is not necessarily true that given equation is integrable even if it 
passes the SC test. (Note that there is no example of integrable difference equation that does not pass the SC 
test.) They have proposed to use the criterion of algebraic entropy instead. Developing the notion of algebraic 
entropy for ultradiscrete equations, together with the development of symbolic manipulation package to deal 
with max-plus algebra, may be quite interesting and important. 
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FIG. 3. Plot of Aio„ to (|l|) for X = 0, X' = p, and Xi = X{ = 1. 
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FIG. 4. An orbit of (jl|) for a = 3, H = 10, X = 11, Xi = 17. Period is 132. 
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FIG. 5. A logarithmic plot of a 132n to @ for a = 3, # = 10, X = 11, X' = 11 + p, Xi = X( = 17. Since slope is about 
25, a 132 „ ~ e 25 ". 
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